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ABSTRACT 

We examine possible phenomenological constraints for the joint evolution of supermassive 
black holes (SMBH) and their host spheroids. We compare all the available observational 
data on the redshift evolution of the total stellar mass and star formation rate density in the 
Universe with the mass and accretion rate density evolution of supermassive black holes, 
estimated from the hard X-ray selected luminosity function of quasars and active galactic 
nuclei (AGN) for a given radiative efficiency, e. We assume that the ratio of the stellar mass 
in spheroids to the black hole mass density evolves as (1 + z)~ a , while the ratio of the stellar 
mass in disks + irregulars to that in spheroids evolves as (1 + z)^ 13 , and we derive constraints 
on a, (3 and e. We find that a > at the more than 4-sigma level, implying a larger black hole 
mass at higher redshift for a given spheroid stellar mass. The favored values for (3 are typically 
negative, suggesting that the fraction of stellar mass in spheroids decreases with increasing 
redshift. This is consistent with recent determinations that show that the mass density at high 
redshift is dominated by galaxies with irregular morphology. In agreement with earlier work, 
we constrain e to be between 0.04 and 0. 1 1, depending on the exact value of the local SMBH 
mass density, but almost independently of a and (3. 

Key words: black hole physics - galaxies: active - galaxies: evolution - galaxies: nuclei - 
galaxies: stellar content - quasars: general - cosmology: miscellaneous 
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1 INTRODUCTION 

Recent work has shown that supermassive black holes (SMBH) are 
ubiquitous in the centers of nearby galaxies. The observational ev- 
idence indicates that the mass of the central black hole is corre- 
lated with spheroid luminosity and mass (e.g., Kormendy and Rich- 
stone 1995; Magorrian et al. 1998) and also with the velocity dis- 
persion of the spheroid (Ferrarese & Merritt 2000; Gebhardt et al. 
2000; Tremaine et al. 2002), suggesting that the process that leads 
to the formation of galactic spheroids must be intimately linked 
to the growth of the central SMBH. In addition, studies of active 
galactic nuclei (AGN) hosts in the local universe with the Sloan 
Digital Sky Survey (SDSS) have also shown that AGN activity is 
closely related to star formation in local galaxies (Kauffmann et 
al. 2003; Heckman et al. 2004). From a theoretical perspective, 
several groups have attempted to investigate the link between the 
cosmological evolution of QSOs and formation history of galax- 
ies within the context of semi-analytic and numerical models of 
galaxy formation and evolution (e.g., Monaco et al. 2000; Kauff- 
mann & Haehnelt 2000; Wyithe & Loeb 2003; Di Matteo et al. 
2003; Granato et al. 2004; Haiman, Ciotti & Ostriker 2004). 

In this paper we make a phenomenological investigation of 
the link between the growth of black holes and the buildup of their 
spheroidal hosts over time. Specifically, we compare the evolution 
of the black hole mass density, pbh, to that of the observed mass 



density of stars in galaxies, p,, and derive constraints on their co- 
evolution. 

To derive pBH as a function of redshift, we use the recently de- 
termined mass density of local black holes (Yu & Tremaine 2002; 
Aller & Richstone 2002; Marconi et al. 2004; Shankar et al. 2004), 
and evolve it back to z ~ 3 using the observed AGN X-ray lu- 
minosity functions (see also Merloni 2004 and references therein). 
The black hole assembly history gathered in such a way assumes 
that the majority of black hole growth is due to mass accretion and 
hence is a function of the efficiency associated with mass to (ra- 
diative) energy conversion. This provides directly the evolution of 
mass accretion rate density as well, given by pBH = 'I'bh. 

To constrain the build-up of stars in the Universe we use the 
independent measurements of the total star formation rate density 
(\&„) and the total stellar mass density (p*) evolution. Because cur- 
rent observations indicate that the black hole mass in galaxies is re- 
lated to the spheroidal component, we will split up the observed p, 
into two terms. We will hereafter refer to the mass density of stars in 
spheroids and disks + irregulars as p sp h and Pdisk+irr respectively. 
In addition to counting the mass contribution from pure spheroid, 
disk, and irregular galaxies, these terms also include masses in the 
spheroid and disk components of individual galaxies. Estimates of 
the ratio of pdisk+irr to p S ph at z — differ by more than a factor of 
5, ranging from 0.2-0.3 (Salucci & Persic 1999; Fukugita, Hogan 
& Peebles 1998; Fukugita & Peebles 2004) to 1.0-1.2 (Schechter 
& Dressier 1987; Benson et al. 2002). The redshift evolution of 
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this ratio should depend not only on the rate of star formation in 
galaxies of different morphologies at different epochs, but also on 
the role of mergers, interactions, and secular instabilities in deter- 
mining galaxy morphology. Both from the theoretical and observa- 
tional point of view, such an evolution has been so far very hard to 
quantify (see e.g. Baugh, Cole & Frenk 1996; van den Bosch 1998; 
Brinchmann & Ellis 2000; Conselice et al. 2004; Combes 2004). 

The rest of the paper is organized as follows. In Section 2 we 
briefly summarize and describe recent measurements of p, and ty* 
as functions of redshift. In Section 3 we derive the black hole mass 
density evolution using the X-ray luminosity function and the lo- 
cal black hole mass function. Under the assumption that p sp h(z) oc 
Pbh(z)(1+z)~ c * and that, at the same time, pdisk+irr/Psph evolves 
as (1 + z)^ 13 , we then compare the history of black hole and total 
stellar mass assembly in the Universe and derive our global con- 
straints on e, a, and j3 (Section 4). Finally, in Section 5 we summa- 
rize and discuss our results and explore their implications. 

Throughout this paper, we adopt a background cosmological 
model in accordance with the Wilkinson Microwave Anisotropy 
Probe (WMAP) experiment. The model has zero spatial curvature, 
a cosmological constant, J1a = 0.71 a Hubble constant Ho = 72 
km s -1 , dominated by cold dark matter with Q m = 0.29 and f2{, = 
0.047 (Spergel et al. 2003). 



2 THE STELLAR MASS DENSITY EVOLUTION 

In this work we make use of all available measurements of the stel- 
lar mass and star formation rate densities from the literature. The 
stellar mass density, p», has been measured in the redshift range 
from z = to z ~ 3 from surveys that select galaxies by rest- 
frame optical or Near Infrared (NIR) light and, in general, estimate 
the stellar masses and mass-to-light ratios of individual galaxies 
by fitting spectral synthesis models to optical-NIR spectral energy 
distributions (see e.g. Brinchmann & Ellis 2000; Cole et al. 2001; 
Dickinson et al. 2003; Fontana et al. 2003; Drory et al. 2004). One 
exception is Rudnick et al. (2003) who have determined the evo- 
lution in p* from the change in the luminosity density and the 
global mass-to-light ratio, the latter being determined by using sim- 
ple models to interpret the cosmically averaged rest-frame optical 
colors. Although cosmic variance is a dominant source of error for 
all but the local values, all support a scenario in which p* at z ~ 3 
was a factor of ~ 7—20 times lower that it is today and that roughly 
50% of the current stellar mass was assembled by z — 1 — 1.5. 

The measurements of p* as a function of redshift also match 
well with the independently measured integrated global star for- 
mation rate density (SFR). The SFR measurements are usually ob- 
tained from e.g., optical emission line measurements, MIR dust lu- 
minosities, sub-mm observations, or observations of the rest-frame 
ultraviolet light. These typically show a rapid rise of the SFR(z) out 
to z < 1.5 where it peaks, followed by a roughly constant level out 
to at least z ~ 4 (see caption of FigureQfor detailed references). 

For the most part, not enough information exists to empirically 
divide p„ and the SFR into different morphological components. 
However, when comparing these data with pbh (see below), it is 
necessary to take into account the fact that local black hole masses 
only correlate with the properties of their spheroidal hosts. We ac- 
count for this by defining a parameter X(z) as the ratio of the mass 
in disks + irregulars to that in spheroids, in the following way: 

p»(z) = Psph(z) + Pdisk+irr(z) = p s ph(z)[l + X(z)]. (1) 

We then assume that \(z) evolves simply according to X(z) = 



Ao(l + z) , where Ao is the value of the disk to spheroid ratio in 
the local universe 1 . 



3 THE BLACK HOLE MASS DENSITY EVOLUTION 

Under the standard assumption that black holes grow mainly by 
accretion, the cosmic evolution of the SMBH accretion rate and 
its associated mass density can be calculated from the luminosity 
function of AGN: <j>(Lt, i, z) = dN/ dL\, a \, where Lboi = eA/c 2 is 
the bolometric luminosity produced by a SMBH accreting at a rate 
of M with a radiative efficiency e. In practice, the accreting black 
hole population is always selected through observations in specific 
wavebands. Crucial is therefore the knowledge of two factors: the 
completeness of any specific AGN survey, and the bolometric cor- 
rection needed in order to estimate Lboi from the observed lumi- 
nosity in any specific band. 

In recent years, enormous progress has been made in under- 
standing both of the above issues. Here we refer to the discussion 
of Marconi et al. (2004), who carry out a detailed comparison be- 
tween optically selected QSOs from the 2dF survey (Boyle et al. 
2000), soft X-ray (0.5-2 keV) selected AGN (Miyaji, Hasinger and 
Schmidt 2000), and hard X-ray (2-10 keV) selected ones (Ueda et 
al. 2003), and work out the respective bolometric corrections. 

The general picture that emerges is that the 2-10 keV AGN 
luminosity function probes by far the largest fraction of the whole 
AGN population. This is primarily due to the fact that hard X-rays 
are less affected by obscuration than either soft X-rays or optical 
light (Marconi et al. 2004). It has also been shown that most of 
the obscured sources, which are the major contributor to the X-ray 
background light, (see e.g. Hasinger 2003, and references therein; 
Fiore 2003) have typically lower luminosity and lie at a lower red- 
shift than the the optically selected QSOs. Therefore, the growth of 
SMBH (and of their associated spheroidal hosts) at z < 1 is best es- 
timated by studying the evolution of the hard X-ray emitting AGN. 
At higher redshifts, optical and hard X-ray selection give consistent 
results (see e.g. Hasinger 2003). 

In the following we will indeed assume that the absorption 
corrected hard X-ray luminosity function of AGN best describes 
the evolution of the entire accreting black holes population between 
z — and z ~ 3. For the sake of simplicity, we will adopt for our 
calculations the luminosity dependent density evolution (LDDE) 
parameterization of the (intrinsic) 2-10 keV luminosity function de- 
scribed by Ueda et al. (2003). 

The redshift evolution of the SMBH accretion rate density can 
then be easily calculated as follows: 

*bh( 2 ) = (1 - £) ^° l(Lx W , z)dL x , (2) 
Jo ec 

where Lx is the X-ray luminosity in the rest-frame 2-10 keV band, 
and the bolometric correction function Lbol(Lx) is given by eq. 
(21) of Marconi et al. (2004), that also takes into account the ob- 
served dependence of the optical-to-X-ray ratio a ox on luminos- 
ity (Vignali, Brandt & Schneider 2003). This in turn yields for the 
SMBH mass density: 



1 Since the stellar mass budget at low redshift is overwhelmingly domi- 
nated by spheroids and disks, Pdisk+irr( 2 = 0) ~ PdiskC- 2 = 0)- At high 
redshift, where irregulars contribute more to the stellar mass density, e.g. 
Brinchmann & Ellis (2000), p d isk+irr ~ Pdisk + Pirregular- 
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pbh(£) = 1 _ r*BH£) dt dz , (3) 

Pbh,o J Pbh,o dz' 

For any given 4>(Lx, z) and bolometric correction, the exact 
shapes of pbh(z) and 'I'bh then depend only on two numbers: the 
local black holes mass density pbh.o and the (average) radiative 
efficiency e. 



4 THE JOINT EVOLUTION OF THE STELLAR AND 
BLACK HOLE MASS DENSITIES 

Given the derived pbh(z) above, we now want to investigate how 
the buildup of spheroidal mass in galaxies is related to the build up 
of black hole mass over time. We make the simple hypothesis that 
the redshift evolution of the spheroid and SMBH mass densities are 
related by the expression: 

Psph(z) = Aopbh(z)(1 + z)~ a , (4) 
where Ao is defined by: 

Ao = p sp h,o/pBH,o = P*,o/[pbh,o(1 + Ao)]. (5) 

For all reasonable values of p*,o, Pbh.o, and Ao, Ao is consistent 
with the expectation from the Mbh — M sp h relation in the local 
universe (Merritt & Ferrarese 2001 ; Mc Lure & Dunlop 2002; Mar- 
coni & Hunt 2003; Haring & Rix 2004). Then, from eqs. (TJ and 
^4) we obtain the desired relation between total stellar and black 
holes mass densities: 

p,{z) = A)PbhO)(1 + z)- a [l + A (l + z)' 13 ]. (6) 

Finally, the total star formation rate ^» as a function of redshift is 
given by: 

».(»)/* = *.w - f^') dxm i ~ z)] 

where x[At(z' — z)] is the fractional mass loss that a simple stel- 
lar population experiences after a time At (corresponding to the 
redshift interval (z' — z)) 2 . For simplicity, we fix z% = 3.5, im- 
plying that the stellar mass density in place then was formed in- 
stantaneously at that redshift. Using the year 2000 version of the 
stellar population synthesis models of Bruzual & Chariot (1993), 
we have found a very good approximation for the function x (f° r 
At> lO 6 ' 5 yr) in the form dx/dt = (At)" 117 . With such a pre- 
scription, the fractional mass loss after 13 Gyr amounts to about 
30%. 

We perform a simultaneous fit to the observed stellar mass 
density and SFR points, assuming the functional forms of J6j and 
above and calculating the black hole mass density evolution ac- 
cording to eqs. {3J and J2j. We vary a, (3, and e over the range 
[0, 1.5], [—2, 1], and [0.02, 0.3] respectively. In FigureQwe show 
the observational data for both p„ (z) and \&» (z), together with the 
best fit curves, in this case calculated assuming Ao = 0.3 and 
Pbh.o = 4.2 x 10 5 M© Mpc~ 3 (Marconi et al. 2004). The shaded 
areas represent the 1-sigma range of models allowed from the joint 
fitting of both datasets. It is important to note that, for different pos- 
sible choices of the local values of pbh.o and Ao (see below), our 
best fits for p* and *D* are almost unchanged and always statisti- 
cally acceptable, with x 2 values typically ranging from 43 to 45, for 

2 An analogous term for pbh. due to the ejection of SMBHs from galaxy 
halos after a merger event, is much more difficult to estimate (see e.g. 
Volonteri, Haardt & Madau 2003) and is neglected here. 
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Figure 1. The top panel shows the evolution in the best fit stellar mass den- 
sity (solid line) as a function of redshift, where the density is given as a 
ratio to the local value, p» j0 = 5.6 X 10 s M Mpc -3 (Cole et al. 2001). 
A value of A = 0.3 and Pbh.o = 4.2 X 10 5 M Q Mpc -3 (Marconi et 
al. 2004) are adopted here. The darker shaded area represents the 1-sigma 
confidence interval for p„. Also shown are the relative decomposition of 
the total stellar mass density into p sp h (dashed line, grey shaded area) and 
Pdisk+irr (dotted line, light grey shaded area), with their corresponding 1- 
sigma confidence interval. Please note that the absolute normalization - and 
to a lesser degree the shape - of these two curves depends on the specific 
value of Ao adopted (see text for details). The data points correspond to 
measurements from the 2dFGRS+2MASS (Cole et al. 2001 - filled circle), 
the Canada France Redshift Survey (Brinchmann & Ellis 2000 - filled tri- 
angles), the MUNICS survey (Drory et al. 2004 - open circles), the Hubble 
Deep Field North (HDF-N; Dickinson et al. 2003 - filled squares), and the 
HDF-S (Fontana et al. 2003 - open squares; Rudnick et al. 2003 - filled 
stars). The lower panel shows the evolution in the best fit SFR (solid line 
and dark red shaded area) and the corresponding black hole accretion rate 
density (rescaled by a factor »4o(l + Ao)). The data points correspond to 
the measurements from Haarsma et al. (2000) - filled black circles; Condon 
et al. (2002) - filled grey circles; Pascual et al. (2001) - open black cir- 
cles; Tresse and Maddox (1998) - open grey circles; Gallego et al. (1995), 
updated following Glazebrook (1999) - filled black squares; Sullivan et al. 
(2000) - filled grey squares; Serjeant et al. (2002) - open black squares; 
Glazebrook et al. (2003) - open grey squares; Brinchmann et al. (2004) - 
filled black triangles; Gallego et al. (1995) - filled grey triangles; Treyer et 
al. (1998) - open black triangles; Gronwall (1998) - open grey triangles; 
Lilly et al. (1996) - filled black stars; Flores et al. (1999) - filled grey stars; 
Cowie et al. (1996) - open black stars; Connolly et al. (1997) - open grey 
stars; Madau, Pozzetti and Dickinson (1998) - filled black diamonds; Pas- 
carelle, Lanzetta & Fernandez-Soto (1998) - filled grey diamonds; Hughues 
et al. (1998) - open black diamonds; Steidel et al. (1999) - open grey dia- 
monds; Sawicki et al. (1997) - black crosses. 
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55 degrees of freedom. Nonetheless, it is apparent from Fig.Qthat 
the best fit solution for p* systematically overpredicts the observed 
points in the redshift range 1.5<z<2, where, however, cosmic 
variance plays a large role and the p» measurements are quite un- 
certain. 

Also shown in Fig. Q are the relative contributions to the to- 
tal stellar mass density of p sp h and Pdisk+irr- It should be stressed 
that such a decomposition is strongly dependent on Ao- The relative 
growth in spheroids and disks + irregulars relative to their respec- 
tive local values is in fact well represented by these curves, as it 
mainly depends on the value of f3, but their absolute growth, and 
hence the total mass density in disks + irregulars relative to that in 
spheroids at any z, is highly dependent on the adopted value of Ao. 
For example, if Ao = 1.0, the dashed and dotted lines in Fig.Qwill 
be moved vertically to attain the same value at z = 0. 

In Figure|2|we show the results as probability contours in the 
2-D parameter space (e, a), obtained by marginalizing the total 
probability distribution over /3, assuming a flat prior on its distri- 
bution, to be conservative. The results are shown for two possible 
values of the parameter Ao (0.3 and 1). As it is well known, a con- 
straint on the average radiative efficiency of accreting SMBH can 
be obtained by comparing the total mass of the relic population 
with the integrated light from all AGN at all redshifts (Soltan 1982; 
Yu & Tremaine 2002; Marconi et al. 2004). Given a reliable inven- 
tory of all the light emitted during the accretion induced growth 
phase (here assumed to be given by the HXLF of AGN), the fi- 
nal result will depend crucially on the value of the local SMBH 
mass density: the lower pbh,o is, the higher the radiative efficiency 
has to be in order not to over-produce the locally measured den- 
sity of relic black holes. This is clearly illustrated by Figure [2] 
that shows the calculated confidence limits assuming two differ- 
ent priors for pbh,o- In one case, shown by the rightmost (b) con- 
tours, we have adopted a lower value as derived, for example, by 
Yu & Tremaine (2002) of p BH ,o = 2.5 x 10 5 M Q Mpc" 3 . The 
radiative efficiency is constrained, at the 3-sigma level, between 
0.08 < e < 0.11. On the other hand, if the local SMBH mass den- 
sity is as high as pbh,o = 4.2 x lO 5 M0Mpc -3 , as suggested by 
Marconi et al. (2004), the radiative efficiency is allowed to take 
a lower value and is constrained, at the 3-sigma level between 
0.05 < e < 0.07, as shown by the leftmost (a) contours. In fact, the 
range of allowed radiative efficiencies would have been larger than 
those shown in Fig.|2|had we also considered the uncertainties on 
the measures of the local black hole mass density 3 . It is clear that 
there is almost complete degeneracy between e and pbh.o, once the 
luminosity function of AGN is fixed. Efficiencies lower than 0.04 
are excluded at the more than 4-sigma level, whatever the prior on 
the local SMBH mass density, provided that it is not higher than 
6.5 x 10 5 MqMpc -3 . This is easily understood, as black holes ac- 
creting with e < 0.04 build up, in the redshift interval between 
z = and z — 3, a mass larger than this value. 

Figure [3] shows the probability contours in the 2-D parame- 
ter space (/3,a), obtained by marginalizing the total probability 
distribution over e, assuming a flat prior on its distribution, once 
again calculated for two different values of Ao. We note here that 
the constraints on these two parameters are almost independent of 



3 We note here that the main difference between the two values of Pbh.o 
comes from a higher normalization of the M-cr relation by a factor 1.6 
adopted by Marconi et al. (2004). It is beyond the scope of this letter to 
discuss thoroughly the reason of this discrepancy. The reader is referred to 
Yu and Tremaine (2002) and to Marconi et al. (2004) for such a discussion. 
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Figure 2. The 1 , 2 and 3 sigma confidence levels for the parameters a and e 
obtained by simultaneously fitting the redshift evolution of the total stellar 
mass density and of the total star formation rate density with eqns. JS} and 
171 . respectively, marginalizing over the other parameter, f3, assuming a flat 
prior on its probability distribution between — 2 < /3 < 1. Four sets of con- 
tours are plotted. Solid lines correspond to a local ratio of the disk + irreg- 
ular to spheroid stellar mass density of Ao = 0.3; dashed to Ao = 1.0. Left 
(a) contours are for pbh.o = 4.2 X 10 5 M Q Mpc _3 (Marconi et al. 2004), 
while the right (b) ones correspond to Pbh,0 = 2.5 X 10 5 MqMpc -3 (Yu 
and Tremaine 2002). 



the adopted value for pbh.o- The two sets of contours in figure [3] 
can be understood as follows. At low redshifts, >]/bh <C 1, and 
Pbh ~ const. Therefore, from eq. {6j, we see that the observed 
decline of the total stellar mass density, and the corresponding in- 
crease of the star formation rate density, will be driven by the term 
(1 + z)~ a [l + Ao(l + z) ]■ For low values of z, this term can 
be approximated by (1 + z) -7 , with 7 ~ a + /3\o/(l + Ao), 
and fitting the low redshift points (see below) with this simple ex- 
pression we obtain the relation a + /3Ao/(l + Ao) ~ 0.3, which 
indeed determines the approximate direction along which the best 
fit parameters a and (3 are allowed to move. 

As it is expected, the lower the value of Ao, the less con- 
strained will the parameter f3 be. A constant \(z) {J3 — 0) is 
allowed at the 3-sigma level only if Ao < 0.3, but the best fit 
value always correspond to a negative f3, i.e. to an increase of the 
disk/irregulars stellar mass fraction with respect to the spheroid one 
with increasing redshift. 

To further test the sensitivity of our constraints on a, (3 and 
e we have also performed a fit restricted to the redshift interval 
< z < 1, where the observational data points are more numerous 
and have smaller error bars. In doing this, we have found that while 
the constraints on a and f3 remain virtually unchanged, those on the 
radiative efficiency are much looser. In particular, for high values of 
the local SMBH mass density, we find that 0.03 < e < 0.12 at the 
3-sigma level. This can be understood because the constraint on e 
comes primarily from the epoch in which the black hole accretion 
rate density is highest, which occurs at z > 1, where the QSO 
luminosity function peaks. The data at z < 1 don't serve to strongly 
constrain e. 
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Figure 3. The contours show the 1 , 2 and 3 sigma confidence levels for the 
parameters a and f3 obtained by simultaneously fitting the redshift evolu- 
tion of the stellar mass density and of the star formation rate density with 
eqns. |5| and 0, respectively, and then marginalizing over the radiative ef- 
ficiency of black holes, e, assuming a flat prior on its probability distribution 
between 0.02 < e < 0.4. Solid contours correspond to the case in which 
we fix the local ratio of disk + irregular to spheroid stellar mass density to 
Ao = 0.3, dashed contours correspond the the case of Ao = 1. The two 
parameters are degenerate along the line a + Ao/3/(l + Ao) ~ 0.3 (see 
text for details). 



In summary, irrespective of all the uncertainties in the determi- 
nation of Ao and in the high redshift stellar mass and star formation 
rate points, the data do indeed put some constraints on the redshift 
evolution of p sp h/ Pbh- In particular, a constant ratio, a = is ex- 
cluded at a more than 4-sigma confidence level, implying that the 
black hole to spheroid mass density ratio must have been larger in 
the past. 



5 DISCUSSION 

It is well known that both AGN activity and SFR decline at low 
redshift. Moreover, the correlation observed locally between su- 
permassive black hole masses and global properties of their host 
spheroids are suggestive of a fundamental link between the forma- 
tion and evolution of black holes and galaxies. 

In this letter we have attempted to constrain phenomenolog- 
ically the joint evolution of supermassive black holes, their host 
spheroids and the total stellar mass density in the universe. 

We have assumed that the black hole mass density evolution 
is due to accretion only and that it can therefore be reliably recon- 
structed from the local black hole mass density and from the red- 
shift evolution of the hard X-ray luminosity function of AGN/QSO, 
with the (average) radiative efficiency of accretion, e as the only 
free parameter. We have then examined the relation between pbh 
and p S ph, assuming that the ratio of these two densities simply 
evolves as (1 + z)~ a . Finally, we have considered the simple 
case in which the ratio of the stellar mass density in disk com- 
ponents + irregular galaxies to that of the spheroid component of 
galaxies evolves as \(z) oc (1 + z)~^ '. By fitting simultaneously 



all the available observational data on the total stellar mass density 
and star formation rate density as a function of redshift, we have ob- 
tained a set of constrains on the three model parameters (e, a, (3), 
which themselves depend on the local values of pbh and A. 

We have shown that the constraint we obtain on the radiative 
efficiency of accretion crucially depends on the adopted value for 
the local black hole mass density, but not on a or /3. 

More importantly, we have shown that BH accretion does not 
exactly track either the spheroid nor the total star assembly (i.e. 
both a and 7 = a + Ao/3/(l + Ao) are larger than zero). This is the 
most important result of our study: irrespective of the exact mass 
budget in spheroids and disks + irregulars, the ratio of the total 
or spheroid stellar to black hole mass density was lower at higher 
redshift. The simultaneous fit we obtain for both p* and ^*also 
indicates that the black hole growth rate is suppressed with respect 
to the total SFR at z < 2, as it can be seen in Figure Q (as also 
shown in the cosmological simulations of Di Matteo et al. 2003). 

Our findings have a direct observable consequences for the 
expected redshift evolution of the Mbh — M sp h relationship: a > 
implies a larger black hole mass for a given host spheroid mass at 
higher redshift. 

We can also compare our results with recent semi-analytical 
work that tries to incorporate SMBH growth and feedback into a 
galaxy evolution scenario. The simplest way to relate black hole 
mass and global galactic properties is via an energy argument, as 
discussed, for example, in Wyithe and Loeb (2003): a SMBH will 
stop growing as soon as the energy provided by accretion onto it 
at a fraction r\ of the Eddington rate for a time tq will equal the 
binding energy of the gas in the DM halo: MBHVLEdd,otQFQ = 
iflbf2mAfhaioUc ■ Here, the black hole mass is in solar units, 
^Edd,© is the Eddington rate of a one solar mass black hole, Fq 
is a coupling coefficient that can be determined from the normal- 
ization of the Mbh — o relation, f2 m and f2b are the matter and 
baryon densities, respectively, and M^aio v c are the mass and cir- 
cular velocity of a DM halo. Most of the feedback models in the 
recent literature (Kauffmann and Haenelt 2000; Cavaliere & Vit- 
torini 2002; Menci et al. 2003; Wyithe and Loeb 2003; Granato et 
al. 2004) assume (sometimes implicitly) a linear proportionality be- 
tween the quasar lifetime and the dynamical time of the dark matter 
halo in which it is embedded: £q = id yn oc £j(z) _1 ^ 2 (l + z)~ 3//2 , 
where £ (z) is a weak function of z and depends on the cosmolog- 
ical parameters only (see e.g.; Barkana and Loeb 2001 for more 
details). Indeed, Wyithe and Loeb (2003) have shown that, if the 
QSO lifetime is equal to the dynamical time of a DM halo, then the 
Mbh — o relation does not change with redshift, while the ratio of 
the black hole to (spheroid) stellar mass density should evolve as 
Psph/PBH oc [£(z)]- 1/2 {l + z)~ 3/2 ~ (1 + Z )- 11S , where the 
final approximation is valid for < z < 2. Comparing this with 
the constraints on a we have obtained in section [4] we see that, 
although marginally consistent with these simple semi-analytical 
schemes, and only for the high Ao case, the data clearly suggest 
a slower evolution for the spheroid to black hole mass ratio than 
that predicted by just assuming a constant Mbh — o relation at 
all 2. Therefore, this would imply that the normalization of the 
Mbh — a relation should evolve with redshift. A caveat is how- 
ever in place here, as of course r\ and Fq in prescriptions such as 
those of Wyithe and Loeb (2003) may not be constant, but vary as 
a function of black hole mass, or mass accretion and as function of 
redshift (e.g.; Merloni 2004). Observational tests of the Mbh — a 
and Mbh — Af sp h relations at higher redshift (see e.g. Shields et al. 
2003) will be crucial to our understanding of the interplay between 
black holes and their host galaxies. 
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Finally, our results suggest that the fraction of p* locked up 
into the non- spheroidal components of galaxies and in irregular 
galaxies should increase with increasing redshift. This result is con- 
sistent with Brinchmann & Ellis (2001) that show that the mass 
density in irregulars increases rapidly out to z ~ 1. 

If we relax the assumption that black hole mass and spheroid 
mass densities are related at all redshifts, then p sp h and pdisk+irr 
in eq. Q would simple correspond to the stellar mass that ends 
up in spheroids and disks, respectively, at z — 0, where we know 
that black hole mass and spheroid mass are related. Under this as- 
sumption, negative values of (5 would imply that the stars in lo- 
cal spheroids were formed at a later time than the stars in present 
day disks. This would directly contradict the relative observed ages 
of disks and spheroids. On the other hand, no such contradiction 
occurs when using our adopted assumption, i.e. that there is a 
Mbh — M sp h relation at all redshifts and that p sp h and pdisk+irr 
correspond, at any given redshift, to the mass in spheroids and 
disks + irregulars respectively. This further supports the idea that 
the dynamical events that lead to the assembly of spheroid stellar 
mass (mergers, secular instabilities, gas accretion, etc.) are directly 
linked to the major accretion events onto the SMBH. 

The upcoming results from projects such as the GEMS 
(Galaxy Evolution from Morphological Studies; see Rix et al. 
2004) survey, that will provide morphologies and structural pa- 
rameters for nearly 10,000 galaxies at 0.2 < 2 < 1.1 and their re- 
spective contributions to the mass density, will surely improve our 
understanding of this issue, and will possibly help us to better un- 
derstand how the simultaneous evolution of black holes, spheroids, 
and disks in galaxies proceeded. 
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